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CN \ We investigate whether it is possible to store and retrieve the intense probe 

O-i' pulse from a A-type homogeneous medium of cold atoms. Through numerical 

^ ■ simulations we show that it is possible to store and retrieve the probe pulse 
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which are not necessarily weak. As the intensity of the probe pulse increases, 
the retrieved pulse remains a replica of the original pulse, however there is 
overall broadening and loss of the intensity. These effects can be understood 
^ ' in terms of the dependence of absorption on the intensity of the probe. We 

On ■ include the dynamics of the control field, which becomes especially important 

. as the intensity of the probe pulse increases. We use the theory of adiabatons 

0\ \ [Grobe et al. Phys. Rev. Lett. 73, 3183 (1994)] to understand the storage 

\ and retrieval of light pulses at moderate powers. 
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I. INTRODUCTION 



Photons are ideal carriers of information. However it is extremely difficult to store 
photons for a long time and to retrieve them from a medium. The basic principle of photon 
storage is based on ultraslow group velocity of the hght [1-5]. The later is made possible 
by electromagnetically induced transparency(ElT) [6] where an external control field is used 
to make an atomic medium transparent even for a frequency near atomic resonance. Under 
such conditions, a probe pulse at a particular frequency and polarization can propagate with 
a substantially reduced group velocity. Experimental and theoretical studies [7-9] show that 
the ultraslow group velocity of the pulse causes spatial compression of a probe pulse that 
has a spatial length of several kilometers in free space. This spatial compression leads to 
localization of the pulses within the atomic medium. Within the spatially localized pulse 
region, atoms are in a state of superposition characterized by the amplitudes and phases 
of the control and probe fields. Therefore, the induced transparency is associated with a 
substantial reduction of the group velocity of the probe pulse which is formed due to a 
coupled field -spin excitation, called dark state polaritons [9]. Switching off of the external 
control field converts the dark state polaritons into purely atomic coherence which is confined 
within the medium. By switching on the external control field at a later time, the atomic 
coherence can be transferred back into the radiation field. The regenerated radiation field 
can be an exact rephca of the original one [10,11] . Therefore, the signal pulse information 
can be stored and retrieved by switching off and on the external optical field [7,9-12]. The 
mode of switching can be adiabatic [9] as well as nonadiabatic [12]. The past studies [7,9-12] 
on this problem have been carried out in the linear regime, where the probe pulse is much 
weaker than control field. Now the question is, whether it is possible to store and retrieve 
an intense probe pulse, in which case, the non-linearity of the medium with respect to the 
probe pulse amplitude becomes important? In this paper we address this question. We 
characterize the probe pulse propagation for different types of pulses in a homogeneously 
broadened A- type medium. The paper is organized as follows: In Sec. 11, we obtain the 
Maxwell-Bloch equations that governs the propagation dynamics of the optical pulse at 
moderate power inside the A-type homogeneous medium. We include the radiative decay 
in simulations. The results of our numerical simulations which also includes the dynamical 
evolution of the control field are presented in Sec. 111. In Sec. IV, we show how the results 
of numerical simulations can be understood in terms of the adiabaton theory of Grobe et al. 
[16]. 



II. DYNAMICAL EQUATIONS FOR PULSES PROPAGATION AT MODERATE 

POWERS 

We consider an atomic system where relevant atomic transitions are taken in the con- 
figuration as shown in the Fig. 1. wherein, one of the transition is coupled to a external 
control field and the other transition is coupled to a probe field Ep. We define all fields 
as 

E(z,t) =£(z,t)e-'('"*-^^) + c.c. , (1) 
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where £^ is a slowly varying envelop, u is the carrier frequency and k is the wave number of 
the field. Here we assume the frequencies of the carrier waves coincide with the frequencies 
of the corresponding atomic transitions = uj-^^ and uj,^ = uo^^. 

We are using the Maxwell-Bloch Equation to describe the propagation dynamics of light 
pulses through the atomic vapour. In the slowly varying envelope approximation the tem- 
poral and spatial evolution of the field envelops is governed by 

dg dg 
dG dG . 

— * — * — * — * 

where 2g — 2diz ■ ^p/^ and 2G — 2dx2 ■ ^c/^ a-re the Rabi frequencies of the probe and drive 

— * — * 

fields, di2, and di2 are the dipole moments, p^g and p^^ are the density matrix elements of the 
corresponding atomic transitions, and c is the velocity of light in vacuum. The propagation 
constant is taken to be equal for all transitions and is given by r/ = SX^N'-f /Sn. It depends 
on the density of the atomic vapor N and the atomic transition wavelength A. 

Within the rotating wave approximation the interaction between the atoms and the fields 
can be describe by the Bloch equations in a rotating frame 

Pii = -2(71 + 72)Pn + iGp^, + igp,^ - iG*p^^ - ig*p^^ , 
P22 = 272p,i + tG*p^^ - iGp^^ , 

P12 = -[71 + 72]Pi2 + iGp^^ + igp^^ - iGp,, , (3) 
Pi3 = -[71 + 72]Pi3 + iGp^^ + i^p33 - igp^^, 

P23 = ^<^Vl3 - WP2I ■ 

Here 7's govern the radiative decay of the state we further assume the equality 71 = 
72 = 7/2. These Bloch equations are to be supplemented by the population conservation 
law 

Pii + P22 + P33 = 1 (4) 

The solution of Eqs.(2)-(3) gives the complete evolution of the atom-field system at any 
instant of time. The analytical solution of the Maxwell-Bloch equations not known though 
under special conditions some solutions are known [15]. Therefore we study the pulse- 
propagation problem only numerically. 



III. NUMERICAL SIMULATIONS 

A. Pulse Propagation: Effect of nonlinearities 

We solve the propagation problem numerically for a homogeneously broadened gas of 
cold A-atoms in the travelling window frame of reference: r — t — z/c, C — We consider 
initial probe pulses of two different shapes given by 
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Gaussian pulse 
\ ^°[secli(^) + / X secli(^)] Sech pulse. 

Here, (7° is the real constant characterizing the peak amplitude of the Rabi frequency before 
the pulse enter the homogeneous medium, a is the band width of the input pulse and Tj 
determines the number of peaks as well as their position. The initial condition for the atomic 
system is to be taken as PsaivC / 1 ^ ^) = 1) with all other density matrix elements equal to 
zero. In order to appreciate the effect of nonlinearities, we first consider the control field 
as a continuous wave G{0,^fr) = constant = 10]. We work under the condition 

of electromagnetically induced transparency. Figures 2(a), 2(b), 2(c) and 2(d) display the 
propagation of two types of probe pulses inside the medium. Prom Figs. 2(a) and 2(c), we 
see that the weaker probe pulse would propagate without any significant absorption and 
broadening inside the medium. The Figs. 2(b) and 2(d) show that the intense probe pulse 
suffers absorption and broadening. The shape of the pulse remains almost identical to the 
input pulse. This behavior of the intense probe pulse can be explained with help of steady 
state probe absorption spectra. In Fig. 3, we show the behavior of the probe absorption 
as a function of the probe detuning when the control field is on resonance. It is clear from 
the Fig. 3 that, increase of the probe field intensity results in the increased absorption of 
the probe for a given frequency in the neighborhood of the frequency satisfying two photon 
resonance condition. Note that the transparency window that appears in the absorption 
spectrum has a finite bandwidth which depends on the intensities of the control and probe 
fields. The width of the transparency window becomes narrow when the probe field intensity 
is increased while keeping the control field intensity constant. Therefore, the condition for 
distortionless pulse propagation is that the bandwidth of the probe pulse should be contained 
within the transparency window of the EIT medium. If the pulse becomes too short, or its 
spectrum too broad relative to the transparency window of the EIT medium, absorption 
and higher order dispersion need to be taken into account. The group velocity dispersion of 
the medium causes the broadening of the probe pulse. 



B. Storage and Retrieval at Moderate powers 

The light pulse propagating in a homogeneously broadened A-type medium, suitably 
driven by another control field, can be stooped and later released in a controlled way 
[7,9,13,14]. The ultra slow group velocity of the fight is the main key issue of the "light- 
storage" technique. In particular, changing the control field intensity results in changing 
the group velocity of the light pulse. The smooth switch off and on of the control field is 
made possible by gradually varying the intensity of the control field with respect to time. 
Therefore, the switching off and on of the control field can be modelled by a super-Gaussian 
shape given by 

Q(Q ^\ _ qO\i e~('^~^)"l ^ ~ ^ adiabatic switching 
^ ^ ^ ' a = 100 nonadiabatic switching. 

In Fig. 4(a), we consider the adiabatic switching of the control field . Switching off of 
the control field give rise to the absorption of the probe pulse when the entire probe pulse 
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is inside the medium. The group velocity of the probe pulse is reduced to zero and its 
propagation is stopped by switching off the control field. The stored probe pulse can be 
retrieved by switching on the control field. The time difference between switching off and 
on is dependent on the life time of the atomic coherence between the state |2) and |3). As 
seen from Figs. 4(b) and 4(d), for weaker probe pulse, the shape of the retrieved pulse is 
same as the original one because the width of the probe piilsc spectrum is very much less 
than width of the EIT window. Therefore, almost perfect storage and retrieval of light is 
possible by adiabatic switching of the control field as pointed out by Fleischhauer et al. 
[9]. When the probe field intensity is large, we observe from Figs. 4(c) and 4(e), that the 
retrieved probe pulse suffers absorption as well as broadening because of narrowing of the 
width of EIT window. Remarkably, the probe pulse can be retrieved even for probe that is 
not necessarily weak. We also present Figure 5 to show the behavior of the atomic coherence 
P32 as a function of retarded time, which depends on the behavior of the input probe pulse. 
In presence of control field, the temporal shape of the atomic coherence is same as the 
shape of the input probe pulse. Remarkably enough the A system leads to storage and 
retrieval even at moderate powers of the probe field. 

C. Nonadiabatic Results 

Scully et al. [12] have shown that for any switching time of the control field, an almost 
perfect storage and retrieval of weak probe pulse is possible. We extend storage and retrieval 
of the probe pulse in the nonlinear regime. The results have been shown in Fig. 6 for an 
intense probe pulse and nonadiabatic switching of the control field. For both adiabatic and 
nonadiabatic switching, the retrieved intense probe pulse is same as original one. However, 
there is overall broadening and loss in intensity of the retrieved probe pulse. The Fig. 6(b) 
shows that the drop in intensity ratio of the output to input probe pulse as function of the 
input pulse intensity. 

D. Dynamical Evolution of the Control field 

The dynamical evolution of the control field becomes important when the intensities of 
the control and probe fields are of comparable strength. Figure 7 depicts time evolutions of 
the control field at different distances. It is evident from this figure that a dip and a bump 
develops in the amplitude of the control field as its propagates through the medium. The 
shape of the bump and dip in the control field depends on the initial shape of the input 
probe pulse at the entry surface of the medium. The propagation dynamics of the dip of the 
control field and broadened probe field together can be understood in terms of adiabaton 
pair [16]. 

IV. ADIABATON THEORY AND ITS RELATION TO LIGHT STORAGE 

In a remarkable paper Grobe et al. [16] discovered what they called as adiabatons. These 
are the pulse pairs which are generated in a A-system under conditions of adiabaticity. We 
show the deep connection of the problem of storage and retrieval of pulses to the theory of 
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adiabatons. The control field is switched on before the probe field, to keep the system in the 
dark state and which is an essential condition for the adiabaton formation. Under conditions 
of negligible damping, the response of the medium then can be very well approximated by 
the solutions 



i d 

i d 
Oio pa — -— 
V dr 

gG 

T/2 



(5) 



P32 



provided the following adiabacity relation is satisfied by the two field: 

dG , . 



(6) 



where V'^ — (G^ + g'^). By inserting solution(5) into the Maxwell equation(2), we obtain a 
pair of coupled nonlinear wave equations 



(7) 



Note that the radiative decay constant does not play any role in the above two equations. 
These two, one dimensional PDE are nonlinearly coupled through the variable V. With the 
help of equation (7), one can easily show that V does not depend on C, implying that V is 
remains constant during the propagation. 



dg 


r] d 




dC 


Vdr 


(f 


dG 


T] d 






Vdr 


iv 



V 



,7T =V{Q,^t) 



(8) 



Thus the conservation law would imply that V any change in probe field is compensated by 
a corresponding in change the control field. Analytical solutions of equation (7) can also be 
obtained by changing the variable r to z{'yT) = \ jT^o ^^(0, ^T)d{'yr): 




y(0,7T)F, 



G{—,1t] =y(0,7T)FG 
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z{-fT) 
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(9) 



-^sW = fi'iO, -2~^(a^)]/^[0, and z^^{x) denotes the inverse function of z. We have 
chosen the initial fields strong enough to ensure the formation of an adiabaton pair. The 
input fields g and G are chosen such that V is constant after a certain time T; hence, for 
T > T the integral z{'yT) can be analytically performed. First we take the control field of 
a constant amplitude (G/7 = 3.16) and the probe field as Gaussian pulse, we then obtain 
(for T>T) 
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9 1 Y'^'^ 



) 






(7r-9o)^) 

— 7^ 



(10) 



2(-fT rTQ- 



V[9''e ^ 



Next the control filed is taken as a super Gaussian shaped pulse and the probe field as Gaus- 
sian pulse as before. The solutions of equation (7) for both these cases, are superimposed 
on the numerical results obtain from density matrix formalism in Figure 7. It is remark- 
able that, the solutions of equations (7) obtain under adiabatic approximation, matches 
extremely well with the numerically solutions of the density matrix formalism containing 
seven nonlinear equations. The adiabatic calculation of atomic coherence P32 of the equa- 
tion(5) is indistinguishable from the density matrix formalism as shown in Fig. 5 [17]. It 
is evident from the temporal profiles of the control and probe at different propagation dis- 
tances that a dip and a bump develops in the control field intensity as it propagates through 
the medium. Figure 7 unambiguously confirms that the adiabaton pair ( consisting of the 
dip in the pump and broadened probe) travels loss-free distances which exceed the weak 
probe absorption lcngth(here typical value of = 2400)by several orders of magnitude 

with an unaltered shape. A gradual decrease of the control field intensity (G/7)^ results in 
a corresponding decrease of the intensity of the probe pulse {9/^Y maintaing the constancy 
of as shown in Fig. 7(b). Therefore, it is clear that when the control field intensity 
becomes zero the probe pulse gets stored inside the medium. The reverse phenomena of the 
retrieval of the probe pulse is achieved by gradual increase of the control field intensity. The 
fact that the numerical results from density matrix formalism matches extremely well with 
the adiabatic approximation clearly indicates that the storage and retrieval of light can be 
understood in terms of the adiabaton pair propagation. 



We have investigated and answered in affirmative, the possibility of storage and retrieval 
of moderately intense probe pulses in a system with A configuration. The propagation of two 
type probe pulses is first analyzed numerically, using Maxwell Bloch equation. It was found 
that the intense probe pulses are absorbed and broadened due to the non-linear dependence 
of susceptibility in the medium. In addition, for the constant control field intensity, the 
width of KIT window becomes narrower with a simultaneous increase of the probe field 
intensity. Therefore, when the intensity of the probe pulse is large, the non-linear effect 
of the medium becomes more important. For larger intensity probe pulses, storage and 
retrieval is possible for both adiabatic as well as nonadiabatic switching of the control field. 
We studied the dynamics of the control field and found that the dip of the control field 
and broadened probe pulse propagate together inside the medium as an adiabaton pair. We 
show that retrieval probe pulse is the part of propagating adiabaton pair. Hence, the storage 
and retrieval of light can be clearly understood in terms of adiabaton pair propagation. The 
same was also studied by an adiabatic approximation. Remarkably, the result of the both 
cases are indistinguishable from each other in the wide parameters range studied here. 



V. CONCLUSIONS 
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APPENDIX: NUMERICAL INTEGRATION PROCEDURE 

The integration of the equation is performed by the fourth order Runge-Kutta method 
[18] .In our numerical calculation, sampling point along r and ( direction are 4 x 10^ and 
2 X 10^ respectively. We prefer the parallezation of the sequential code because of the large 
number of sampling points which requires large execution time. The sequential code can be 
easily parallelized if each iteration is independent of the other, that is no variables that are 
written in some iteration will be read and/or written in another iteration. But in our code 
each iteration is dependent on the previous iteration. Thus for our flow dependence case, 
the parallelization of the code is difficult and we do the parallelization manually by using 
OPEN MP directives in RS6000 in IBM machine. 
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FIG. 2. The probe field intensity in the medium is plotted against retarded time at different 
propagation distances within the medium. Fig (a) and (c) show the probe pulse propagation with 
non-diminishing amplitude, for small intensities. Fig (b) and (d) depict the broadening and loss 
of intensity in case of an intense probe pulse case. In all the cases the control field is taken as 
CW(G = 3.I67). 
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FIG. 3. Imaginary parts of susceptibility [Xis^^Z-^l^^is P] ^ probe frequency in the pres- 
ence of control field as a CW G = 3.1 67. The width of the transparency window decreases with 
increase in the intensity of the probe field. 
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FIG. 4. (a) shows the intensity of the control field as a function of retarded time at the entry 
surface of the medium at ( = 0. Switching mode of the super Gaussian control field is adiabatic. 
The frames (b) and (d) show the time evolution of the weaker probe pulse at different propagation 
distances; and the frames (c) and (e) depict the time profile of the intense probe pulse at different 
propagation distances. 
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FIG. 5. Fig. (a)-(d) shows the temporal profile of atomic coherence — P32 against retarded time 
at different propagation distances. 
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FIG. 6. (a) shows storage and retrieval of intense probe pulse even for nonadiabatic switching of 
the control field, (b)drop in intensity ratio of the probe retrieved to the input pulse as a function 
of input probe intensity for the case of non adiabatic switching of the control field ; the lout is 
measured at t/C/t = 3200 and 7r = 1000. 
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FIG. 7. (a) and (b) shows temporal profiles of the control and probe field (5/7)^ 

at different propagation distances within the medium. In Fig. (a) the input control field 
is a CW. In Fig (b) the input control field is a Super-Gaussian shape with parameter 
T2 = 575/ J, af = 200/7. The common parameter of the above two graphs are chosen as: 
G° = 3.167,5° = 1.147, TO = 200/7,0- = 90/7. The results of simulations using Maxwell-Bloch 
equations are indistinguishable from the results based on adiabaton theory. 
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